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ABSTRACT 


This work simulates an acousto-optic signal processing technique 
that can be used to filter radio-frequency (RF) electronic signals. Of 
particular interest is the removal or excision of narrowband 
interference from broadband signals. The program developed 
simulates the performance of a binary phase shift keying (BPSK) 
signal with a narrowband interference. The width and depth of 
notches that model an excisor-detector array are variable. In 
addition there is a capability to move the notches in order to place a 
notch over the interference for maximum filtering capability. <A 


summary of several possible excisor techniques is presented. 
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DISCLAIMER 


The reader is cautioned that computer programs developed in 
this research may not have been exercised for all cases of interest. 
VVhile every effort has been made, vvithin the time available, to 
ensure that the programs are free of computational and logic errors, 
they cannot be considered validated. Any application of these 


program without additional verification is at the risk of the user. 
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I. INTRODUCTION 


A - BACKGROUND 


Many companies are designing, developing, and field-testing 
Acousto-optics Spectrum Analyzers (AOSA) for communications, 
radar and electronic warfare surveillance. In many such 


applications, acousto-optics is not merely a good choice— it is the 


only choice. [Ref. 1] 
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Figure 1 - Acousto-optic spectrum analyzer 


As an introduction to acousto-optic spectrum analyzers consider 


the system in Figure 1. Coherent light illuminates a one dimensional 


spatial light modulator, in this case a Bragg cell. The light diffracted 
by the cell is collected by the lens and focused onto a photodiode 
array. The light amplitude at the photodetector is the spectrum of 
the information in the Bragg cell. This serves to illustrate the 
simplicity of one optical technique for calculating spectra. [Ref. 2] 

The problem of rejecting relatively narrowband interference 
signals is becoming more important as the bandwidth of digital 
communications systems increases. Acousto-optics is showing 
promise in this area [Ref. 1]. An optical implementation of this 
technique is shown in Figure 2. 

In order to implement an optical excisor we need more elements 
interacting in the AOSA basic architecture than was presented in 
Figure l. As showed in Figure 2, the signal s(t) is applied to a Bragg 
cell (1). The cell is illuminated with coherent light) (2) and the 
spectrum of the signal is formed on the photodetector array (3). The 
light at each photodetector has an amplitude. proportional to the 
amplitude of the corresponding acoustic frequency in the Bragg cell. 
A given acoustic frequency component acts as a moving diffraction 
grating; for this reason the light diffracted by it is Doppler 
frequency-shifted by the frequency of the acoustic signal. The 
acoustic signal component can be reconstructed by down-converting 
the diffracted light to the acoustic frequency. This is achieved by 
adding the optical signal to a reference beam (4). The resulting 
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COIHIUIIICG. UCâ aii On a priotioàciector array. The nhotodetector 


to 


produces an electrical signal vvith a frequency equal to the difference 
between the frequencies of the two optical inputs. This results in the 
recovery of the acoustic frequency components and hence the 
recovery of the input radio frequency component. 

The optical excisor continuously forms the spectrum of a 
broadband electrical signal and reconstructs an exact replica of the 
input signal from that spectrum. The important feature of this 
process is that the spectrum is spatially distributed and is physically 
available at the detector plane for removing an undesired portion of 
the spectrum. For example, the excision of narrowband interference 


from broadband signals is possible. 
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Figure 2 - Optical excisor architecture 


Two significant characteristics of an optical excisor are: 1) it has 


an RF signal input and an RF signal output, 2) it is a predetection 


signal processor producing a continuous coherent output [Ref. 3]. 


Optical Excisor 





Figure 3 - Optical excisor characteristic 


B- OPTICAL EXCISION 


If we have two input signals that are sufficiently different in 
frequency to be resolved in the back focal plane, a spatial light 
modulator can be placed in this location and used to turn off either 
one or both of the input signals [Ref. 3]. 

Another technique would be to block the light of one signal by 
using an aperture or small wire in the focal plane, thereby allowing 
passage of only the desired signal. [Ref. 3] 

A broadband signal that has a narrowband interference signal in 
its spectrum can be applied to the input of the optical excisor. The 
narrowband signal will form a small spot of light, while the 
broadband signal will extend over a relatively large area. Using a 
small wire, the strength of the interference can be significantly 
reduced while only blocking a small portion of the broadband signal. 


The bigger the wire, the more the interference will be blocked, and 


therefore the more it vvill be attenuated at the output. However, a 
large wire will also remove more of the broadband output signal and 
perhaps begin to cause significant degradation in the broadband 
output waveform. The task for the designer of an optical excisor is to 
optimize the transfer function to achieve a useful amount of 
attenuation of the narrowband interference without removing 
significant portions of the broadband signal. [Ref. 3] 

An alternative excision approach is to place a material such as 
photochromic glass or photodichronic crystal in the focal plane. 
These materials will darken in an amount proportional to the 
irradiance of the incident light (much like photographic film, but in 
real time). Narrowband interference signals, because of their 
relatively high intensity, will be attenuated automatically [Ref. 3]. 

As can be seen, using different excision methods in the optical 
Fourier transform plane will create excision notches with different 
properties. 

The following is a summary of the different experiments found 
in references which illustrate various techniques for excision as well 
as their results. 

Thin opaque wires of known diameters placed in the transform 
plane can create narrow excision notches. Using thin bare copper 
wire in the transform plane, it blocks a portion of the diffracted 
signal spectrum, thereby creating an excision notch. In different 


experiments it was found that i enuation achieved was about 40 


dB with inexpensive lenses. [Ref. 6] 


A Pockel's Readout Optical Modulator (PROM) spatial light 
modulator can be used in a fashion which automatically attenuates 
narrowband  interferences. The PROM has nonlinear clipping 
characteristics which attenuates intense signals, while weaker signals 
pass unaffected. In fact, a strong signal can be attenuated below the 
level of a weak signal. Broadband signals should be relatively 
unattenuated compared to more powerful narrowband signals 
because their power spectrum distribution in the optical Fourier 
transform plane will be below the clipping threshold of the PROM. 

The PROM optical clipper is intended to be utilized as an image 
storage device which is written with blue light and read out with red 
light. The experiment reported in reference 3 utilized the red HeNe 
laser light for both writing and reading functions. A strong 
narrowband interference signal will create a bright red spot in the 
transform plane which effectively acts to write a dark spot on the 
PROM. This dark spot prevents the narrowband interference from 
being passed to the photodetector. Thus the strong narrowband 
interference signal effectively attenuates itself. [Ref. 3] 

Reference 3 showed that the PROM-based excisor has attenuated 
a narrowband signal by 28 dB relative to the broadband pulsed 
signal. It demonstrated encouraging results, but the utility in a 
practical system is limited by the fact that a PROM is an image 
Storage device which must be cycled through an erase phase 


followed by a read/write phase. 


Variations of the PROM are being considered vvhich vvill avoid 
the image storage properties of the current PROM so that time- 
framing and cycling problems may be avoided. [Ref. 3] 

A Programmable lead lanthanum zirconate (PLZT) ceramic 
modulator can also be used in the transform plane. A programmable 
PLZT ceramic light modulator will allow the filter magnitude 
response to be controlled by electronic means. By applying a 
transverse electric field to the PLZT ceramic, the polarization of the 
transmitted optical beam can be rotated. By placing the PLZT 
ceramic between crossed optical polarizers, the intensity of the 
transmitted optical beam can be electronically controlled. Spatial 
control of the transmitted optical intensity can be obtained by having 
separately addressable electrodes across the aperture of the PLZT 
ceramic. [Ref. 3] 

The results of the PLZT-based optical excisor showed that 18 dB 
of excision could be obtained while an opaque wire having a similar 
notch width obtained 23 dB of excision depth. [Ref. 3] 

Another form of excision is due to the spaces between detectors 
in a detector array to perform the optical excision. In the frequency 
domain, they are represented by a set of notches that are 
periodically spaced. The width of the notches is proportional to the 
spaces between the detector elements. The depth depends on the 
width and the diffraction effects (A similar model also applies when 
a linear array of fiber optic cables leading to individual detectors is 
operties). [ft is this laiter form of 


excision that we wanted to study. 


C- THESIS OBJECTIVE 


The objective of the research described in this thesis was to 
develop a computer program which, through simulation, would 
provide a capability to evaluate the filtering ability of gaps between 
detectors in a detector array versus narrowband interference in a 
broadband signal. 

Using the Naval Postgraduate School's library and computer 
facilities, we wanted to develop a computer program based in 
MATLAB386 [Ref. 4] using a Compaq 386/20 computer. 

The specific goals were to evaluate the distortion on the original 
signal due to the filtering process and to examine the notch 
dimensions (width and depth) for their effect on the filtering 
capability. 

Chapter 2 serves as a brief description of the optical excisor 
system. This chapter describes an acousto-optic signal processing 
technique that can be used to filter radio-frequency electronic 
signals. A discussion of the mathematical model chosen for the 
simulation is given in Chapter 3. The detailed description of the 
computer program is presented in Chapter 4, showing the 
information available during the simulation. Chapter 5 contains the 
simulation results and an evaluation of the program, and finally 
Chapter 6 contains concluding remarks and suggests areas of future 


research. 


II. GENERAL DESCRIPTION OF THE SYSTEM 


The optical excisor separates the input signal s(t) into each of its 
Fourier components, processes the magnitude and phase of each 
Fourier component and finally sums the modified Fourier 
components into a filtered signal. A general description of an optical 
excisor can be based on the architecture shown in Figure 5. 

The electrical input signal s(t) typically is derived from a 
broadband antenna, downconverter, and amplifiers. This signal is 
applied to a Bragg cell illuminated by a coherent light, illustrated in 


Figure 4. 
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Figure 4 - Bragg diffraction in a spectrum analyzer 


The acousto-optic modulator is a block of transparent material 
with a piezoelectric transducer bonded to one end. The other end is 
tapered to trap the acoustic waves and loaded to provide an effective 
acoustic termination. A piezoelectric transducer is sandwiched 
between electrodes and vibrates in synchronism with the applied 
voltage. When driven by an electrical signal, the transducer excites 
an acoustic which is a replica of that signal. This acoustic wave 
propagates down the long axis of the Bragg cell and produces 
compressions and rarefractions in the material [Ref. 2]. The 
interaction between coherent light and the acoustic wave results 
from an increase in the index of refraction of the material at a 
compression and a decrease in the index at a rarefraction. This 
results in a modulation in phase of the light passing through the 
acousto-optic modulator. The resulting spatial phase distribution can 
be considered equivalent to a superposition of plane waves of 
various phases and amplitudes. By using a lens to focus the various 
diffracted plane waves, each wave can be separated and the 
amplitude and phase of each spectral component of the original input 
signal can be measured. So far the description is similar to an AOSA 
(Ref. 3). 

In addition to diffracting the plane vvaves based on the spectrum 
of the input signal, the acousto-optic modulator also Doppler shifts 
the diffracted plane waves. The amount of Doppler shift is exactly 


equal to the frequency of the input signal, and therefore, if the 
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amount of Doppler shift is measured, the frequency of the input can 
be determined. [Ref. 3] 

An easy way for measuring an optical Doppler shift is to add 
some unshifted light from the optical source and process the summed 
signal electronically. The resulting beat signal frequency will be 
equal to the frequency difference between the two signals and hence 
be equal to the amount of Doppler shift. (Ref. 3] 

The unshifted reference beam, in Figure 5, is derived from the 
spatially coherent optical source using optical splitters (partially 
reflecting mirror). The light reflected by this beam splitter is set to 
mirror 1 which then reflects the light to the acousto-optic modulator. 
The light transmitted by the beam splitter is reflected by mirror 2 
up to a second beam splitter, reflecting some of the reference beam 
to the photodetector array. The second beam splitter also passes 
some of the light diffracted by the acousto-optic modulator. By 
adjusting the mirrors and beam splitters, both the reference beam 
and the diffracted beam appear to come from the same direction and 
are added together by the beam splitter before the detector array. 
By using a spatial light modulator in the transform plane and 
Selectively passing only the desired portions of the spectrum, 
interference signals can be rejected [Ref. 3]. The photomixing and 
signal summing by the detector reconstructs the time domain 
waveform of the desired portion of the input signal and provides a 


coherent continuous output i(t). 


Z 401W rors 


uva y 








ug 32u2342f2 


aunjd 

ui 4ofSuv4 y u! 

40)v9npoJjw 

| 1y3r] joynds 





cessa mm peu 


sindino 7 = (0 


431110 x 
(1)5 puI 






ADIY 
40jJajapoioud 


c" sooo c cece ees 
a 


-—---- 





dois pondo sus] 
I ul40fsuv4 Jj 102]npo w 
E 2ydo-0jsn02y 





SaSUua'] 


duipupdx:] uos] 


I 4041 


oca 





334N0S 14 3Y] 
1u342107) 101 Dd Ç 


Figure 5 - Optical excisor 
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The five main components of the system presented are: 
(1) spatial coherent light source, 
(2) one-dimensional, acousto-optic modulator; 
(3) lens for Fourier transformation; 
(4) a frequency plane spatial modulator, and 
(S) the photodetector array. 

The present work considers the detector array acting as a spatial 
light modulator in the frequency plane; for this reason, from now on, 
it will be called as "excisor-detector" array. This technique is of 
particular importance in rejecting narrowband interference in 
wideband communications systems. By using the spaces between 
detectors as notches in the frequency domain, one can control the 
frequency rejection to match the signal environment. The signal 
environment may consist of a wideband signal with a narrowband 
interference signal located at some arbitrary frequency within the 
bandwidth of the desired wideband signal. The object is to block the 
narrowband interference signal with one of the narrow bandwidth 


freguency rejection notches of the excisor-detector array. 
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III. THE MATHEMATICAL MODEL 


A more precise description of the operation in the optical excisor 
can now be given. The broadband signal including a narrowband 
interfering signal is s(t). The collimated coherent input light is 
diffracted through an angle proportional to the acoustic frequency 
present in the Bragg cell. This acoustic frequency results from the 
applied electrical signal s(t), where s(t) is real. The diffracted optical 
signal amplitude at the exit port of the cell is proportional to s(t-x/v), 
where x is the distance from the input transducer and v the acoustic 
velocity. The acousto-optic (AQ) cell has a finite aperture equal to T, 
where T is the transit time of the acoustic wave in the AO cell. This 
finite aperture corresponds to an apodization of the acoustic signal 
such that the signal may now be represented by s(t-x/v-T/2) [Ref. 2]. 
Also for this analysis, the AO cell is assumed to operate over an 
octave RF bandwidth (or less) and the positive first order diffraction 
is used [Refs. 2, 3 and 5]. 

The optical system takes the Fourier transform with respect to 


the spatial variable x. 


e (1) 


where L is equal to the length of the AO cell and p is the spatial 
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frequency variable. Associated with each spatial frequency p is a 
temporal frequency œ = pv. 


A change of variables (where u=t-x/v-T/2) yields: 


. T 
-jJt——)v rt-T 
S(p,t) =ve "5 1 J s(u) e" “du 
i (2) 


ET ! 
ony = eÌpvu g 
(pv) J s(u) u " 


Equation (2) is seen to be the Fourier transform of the portion of 
s(t) within the aperture of the Bragg cell. The exponential phase 
term indicates that the light diffracted to position p is frequency 
shifted by an amount equal to the temporal frequency o. 

The detectors respond to light intensity rather than light 
amplitude, for this reason it is necessary to bring in a reference wave 
intersecting the frequency plane. The reference wave has a flat 
wavefront which may be parallel with the frequency plane or 
intersect the frequency plane at a small angle ©. This reference wave 
with tilt angle Q may be written as Ael$P, where A is the amplitude 
of the reference wave [Ref. 5]. 

The total amplitude in the output plane, T(p,t), is the sum of the 
reference wave and the transform of the signal multiplied by the 


filter signal, H(p,t), (representing the spaces between detectors). 


| 


The corresponding intensity distribution is: 


2 
I (p.t) = M3) (55 


2 2 2 2 2 
TE Ð| =Hp,t) A +v" S(py) H(p,t) 
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where Re(x) means the real part of (x). The output from the excisor- 
detector array is the summation of all photodiodes; this corresponds 


to an integration over the spatial variable p. 


Í 2 

Lo(p. 9 dp - J [Tip. of dp bj 
el . 2 

- [ luo. t) v E 2) S(pv) 42A e? P Hp, 0| dp 


= Jv? SN H(p, t) dp +A’ SHG, y” dp 


. a 
32A v Re] Hp, 0 S(py) sib * Yt 7 2). qp (8) 


As stated before, the s(t) signal is assumed to be real and to 
occupy an octave or less of bandwidth. Under these assumptions the 
first two terms of equation (8) do not contribute in-band frequency 
components to the reconstruction of the original input signal [Refs. 2 


and 5). 
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The last term of equation (8) represents an inverse Fourier 
transform of the original Fourier transformed signal, which yields the 
original signal with time delay equal to half the time aperture of the 


AO cell plus the delay equal to the reference beam. 
Ie + (t - D) 
2 A vRe | H(p,t) Síp v)e? v 2” dp= 


=2Ah(t)s(t-> + 2) 
(9) 


Equation (9) represents the filtered reconstruction of the original 
signal. 

The mathematical model described can be summarized in Figure 
6 for the program calculations. The test BPSK signal Vi(t) is injected 
in the input. The guantity Vi(f) represents the FFT of the input 
signal, H(f) is the transfer function of the detector-excisor array, and 
Vo(f) represents the filtered signal. Recovering the filtered signal 
back to the time domain, the IFFT is performed to obtain Vir(t) 
which is electronically processed to detect the original binary 


information. 


| 


BPSK (t) Signal 


Vi(t) 
Vi (t) 2 m(t) .cos( 2 Af 1) 
Vi(f)-fftVi(t)) m(t) is a Random Binary 
Voltage [+1 or-1] 
SS FFT 
Vi(f) 
Vo(f)>Vi(f) . H(f) 
V. HO 
Vo(f) 
Vir(t)-ifft( Vo(f)) 
PE IFFT 
Vir(t) 
Vr (t)-Vir(t).cos(2 xft) E Coherent 
i e i K unmodulated 
~~ M carrier 





cos ( 2 aft) 


Recovered Random 
MU Binary Voltage 


Figure 6 - System representation 
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IV. COMPUTER PROGRAM DESCRIPTION 


A - INTRODUCTION 


The purpose of the simulation program is to model and analyze 
the effects of an optical excisor. Of particular interest is the removal 
of narrowband interference from a BPSK communication signal. This 
interference is removed using the notches between detectors, in a 
detector array, as a filter. 

Optical excisor performance is measured in terms of how well 
the offending signal frequency components are removed from the 
output, and how much of the original signal of interest is corrupted 
through loss of frequency components during the filtering process. 
[Ref. 6] 

The excisor-detector array contains gaps between detectors that 
appear as notches in the freguency domain (filter). For our purpose 
the width and depth of the notches can be varied to simulate optical 
and mechanical adjustments and the effects produced on the test 
signal can be observed. 

The modeling is intended to be simple enough as not to obscure 
the fundamental ideas but complete enough to include important 


effects and avoid incorrect conclusions. 


= 
M2 


B- PROGRAM STRUCTURE 


The program was 


divided into eight different blocks 


that 


represent various stages in the process. Each of these blocks 


performs specific functions presented in Table 1. 


A simplified diagram of the program is presented in Figure 7 


which gives the general 


structure of the system. 


TABLE I - BLOCK FUNCTIONS 


FUNCTIONS 


e INP B.m 


e M BPSK.m 


e FILTER.m 


e S BPSR.m 


e FFT_B.m 


e S FIL.m 


e IFFT. B.m 


e COD B.m 


l- Lets you input the data for the simulation. 


1- Master Block executes other sections. 


2- Controls the loops to simulate the movement of 


the Excisor-detector array. 


l- Builds the filter. 


2- Generates the frequency axis. 


l- Creates the BPSK signal. 
2- Creates the interference. 
3- Generates the random numbers. 


4- Generates the time axis. 


l- Calculates the FFT of the input test signal. 


2- Calculates the FFT Theoretical for comparison. 


l- Computing the filtering process in the 


frequency domain. 
l- Calculates the IFFT of the filtered signal. 


l- Coherent detection of the filtered signal. 
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C- THE PROGRAM 


In this section a detailed explanation of each block is presented 


for a better understanding of how the computer program works. 


Filter.m 


Keyboard 





Figure 7 - Program structure 


zal 


1 - BLOCK INP B.m 


The input data block is for introducing the variables needed to 
run the program like the carrier frequency (fc), the amplitude of the 
signal (Ac), the central frequency of the interference signal (fi), the 
step interval in order to simulate the movement of the detector 
array (delta), the number of loops that is needed for the program (I), 
and the seed number for the random number generator (seed). All 


values are presented on the screen. 


2- BLOCK M BPSK.m 


This block controls the sequence during the simulation and the 
loop to simulate the movement of the excisor-detector array. In 
order to reduce the amount of memory, the data is stored in 
temporarily files to be used by the other section. After each block 
has been run and the data has been stored in a temporary file, the 
memory is cleared before executing the next block to keep the 


computer from running out of memory. 


3 - BLOCK FILTER.m 


The main functions of the Filter.m block are to build the filter 


transfer function and to create the frequency axis that will be used 


by the simulation. The filter consists of a signal that represents a 
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detector array in the frequency domain. An ideal representation of 
the excisor-detector array is showed in Figure 8. In Figure 9, the 
filter magnitude response of the model of the excisor-detector array 
is presented. The width (w) may vary from O to 120 KHz and the 
depth (d) from O to -90 dB. Figure 8 also shows the center frequency 
for each notch in the excisor-detector array assuming that the array 
is centered at 70 MHz. Appendix C provides useful information to 
work with the filter, like the freguency at which any notch starts and 
ends for a given width and the correct value of notch depth that 
represents the desired attenuation. 

Performing this RF filtering process in the freguency plane of an 
optical spectrum analyzer has advantages over the electronic 
techniques [Ref. 6]. When an object is placed in the frequency plane 
of the optical spectrum analyzer, the magnitude of the affected 
frequency interval is attenuated and the phase remains constant. 
Also, many filters may be realized simply by varying the size of the 
notch in the frequency plane. 

The excisor-detector is placed at one focal length away from the 
Fourier transform lens into the frequency plane. It has 13 
photodiodes elements that represent the different channels 750 KHz 
wide, adding to a total 9.75 MHz window as shown in Figure 8. 


(These values represent values used in an actual laboratory setup). 
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Figure 8 - Excisor-detector array 
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Figure 9 - Filter magnitude response 


In the computer program (Appendix B, p. 66) the notch depth is 
selected to represent the desired attenuation. The separation 
between points in the frequency domain is 5 KHz. To represent a 
notch 60 KHz wide, 12 numbers or different elements are needed. 

The frequency axis is created considering that the sampling 
frequency (fs) must be greater than twice the highest frequency 
present in the signal that is being sampled. That is, the sampling 
frequency must satisfy the inequality fs>2fmax [Ref. 7]. The highest 
frequency in our test signal is fmax=70 MHz, so the sampling 
frequency must be greater than the Nyquist frequency 2fmax=140 
MHz. For the simulation fs=163.84 MHz has been chosen. 
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Using MATLAB386 (Ref. 4) the following expression was used to 
create the frequency axis: 
f = (bin_number -1)*fs/N (10) 
where 
f = frequency axis, 
fs 


N = number of points or vector length, and 


sampling frequency, 


bin number = N/2+1. 


In order to represent the filter signal centered at 70 MHz the 


values used in the previous expression have been N = 32768, 


bin_number-1 = 16384, and fs = 163.84 MHz. 


Two plots are presented on the screen, showing the detector 


array filter signal and the filter magnitude response. 


4- BLOCK S BPSK.m 


This block constructs the input test BPSK signal, the interference 
signal, and the time axis to be used in the program. 


The BPSK signal is represented by [Ref. 7] 


s(t) = Ac cos[@c t+Dp m(t)] (11) 


where m(t) is a bipolar baseband data signal with peak values +1. 
Recalling that cos(x) and sin(x) are even and odd functions of x, the 


BPSK signal can be obtain by [Ref. 7] 
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A elit 


s(t) = (Ac cos Dp) cos Oct - (Ac sin Dp ) m(t) sin œc t 12) 
first term second term 
where the first term is called the pilot carrier term and its level is 
set by the value of the peak deviation AO = Dp. 
For digital angle-modulated signals, the digital modulation index, 


h, is defined by [Ref. 7] 
h= 2 AB /T (13) 


where 2 A0 is the maximum peak-to-peak phase deviation during 
the time that it takes to send one symbol, Ts. For binary signaling, 
the symbol time is equal to the bit time, Ts = Tb. 

The power in the second term (called the data term) needs to be 
maximized. This is accomplished by letting AO = Dp = 90? [Ref. 7], 
which corresponds to a digital modulation index of h = 1. For this 


case, changing sin by cos, the BPSK signal becomes 


s(t) = Ac m(t) cos Oct (14) 

where 
Wc = 2 1 fc, and 
Ac = Amplitude. 

Equation 14 is the algorithm applied in the test signal 
generation. The m(t) signal is a random binary data that represents 
the modulating waveform originated from a digital information 
source and is illustrated in Figure 10. The use of random data should 


be more effective than using any other deterministic method. 


2a 


9. AA 


it 


-2 SS eS ee D'ara Ges rames 


0 0 2 04 0.6 0.8 


AMPLITUDE, volts 





TIME, scc x10-* 


Figure 10 - Binary random data 


The random number generator (RNG) employed must fulfill the 
following requirements [Ref. 8]: 

e Uniform: It must pass several statistical test for U(0,1) uniform 
distribution. 

* Reproducible: To be useful, the RNG should reproduce the 
same string of numbers given the same seed every time; this 
is of particular importance for testing purpose. 

* Minimum Memory: It should be efficient to use and fast. 

* Non-repeating: A number sequence must not repeat itself 
during a simulation run. 

e Statistical Independent: This property is very difficult to 
achieve because of the recursive properties of the generator. 

The mathematical algorithm was based on the linear 


congruential method of the form (Ref. 8): 
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zi = (a*zi-1 + c) MOD m (15) 
where 
a = multiplier, 
c = increment, 
m = modulus, and 
Zo = seed. 
It can be shown that the pseudorandom number generator can 


be constructed using the following expression [Ref. 8] 


zi = (16807 *zi-1) MOD 231.1 (16) 


For the random case the test signal generated by Eq. 14 will 
produce a continuous spectrum as will be presented later. 

The interference signal is also generated in the present block 
based on a sinusoidal signal with an amplitude three times the 
original BPSK signal (in the time domain) and a carrier frequency 
selectable from the input data. As a result of adding those two 
signals, the test BPSK signal plus interference is produced. 

The generation of the time axis is other function of this block. A 
total time of 0.2 milliseconds with intervals of 6.1 nanoseconds has 
been selected. 

Plots presenting the random binary data, BPSK test signal, and 


BPSK test signal plus interference are presented on the screen as 
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shown in Figures 11 and 12. A carrier frequency of 70 MHz and an 
amplitude of 1 volt was used for the BPSK signal and a central 
frequency of 67.4 MHz and 3 volts amplitude for the interference 


signal. 








0.5 


MIL 

AAA AA 

805) TEL IH PEI i i 
ACC 
IL | v | H 


AMPLITUDE, volts 


TIME, sec x10-9 


Figure 11 - BPSK signal and modulated waveform 


3-SBEOGRSBBT Bm 


The FFT b.m block uses the input signal constructed in the 
previous block and calculates the fast Fourier transform. The first 
step in applying the discrete Fourier transform is to choose the 
number of samples N and the sampling interval T, and, secondly, to 
interpret the results correctly. The discrete transform approximates 
the continuous transform rather poorly for the higher frequencies 
[Ref. 7]. To reduce this error it is necessary to decrease the sample 


interval T and/or to increase N. 


3 U 


In our case the relationship used vvas 


T Period -N (17) 
where 


N = 32768 (N = 2? where @ is an integer value), and 
Period = 0.2 millisec. 


So 1 = 6.1 nanosec. 


AMPLITUDE, volts 





TIME,sec x10-8 


Figure 12 - BPSK signal plus interference 


From MATLAB386 (Ref. 4] the FFT(x) is the discrete Fourier 
transform (DFT) of vector x, computed with a radix-2 fast Fourier 
transform algorithm. If the length of x is not an exact power of two 
then it is padded with trailing zeros. This is not our case because we 
have selected N = 32768 = 2!°. 

The function that implements the transform is given by 

N- 1 


X(k+1)= Y xn+1) we 
pos (18) 
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vvhere Wr = e JN) and N - length(X). Note that the series (18) is 


written in an unorthodox way, with the subscript from k+1 instead of 
the usual k because subscripts on Matlab vectors run from 1 to N 
instead of from O to N-1 [Ref. 4]. 

In order to plot the FFT results, the frequency axis generated in 
the Filter.m block is used. Actually the power spectral density (PSD) 
is calculated and plotted as shown in Figures 13. 

Also, for comparison purposes, the theoretical power spectral 
density of a BPSK signal with the same characteristics as the test 


signal was calculated and plotted using the following algorithm: 


2 


2 NP 
A sitio 
P(f)- TET 


4R n(f—f.)T, 
(19) 


where 
Tp = bit time = 1/R, and 


R = bit rate. 


The resulting BPSK spectrum is shown in Figure 14. The null-to- 
null bandwidth is 2R. The figure shows plots of the calculated and 
the theoretical PSD of the test signal. As we can see, the calculated 
spectrum (solid line) appears to be noisy but follows the theoretical 
spectrum plotted using equation 19 (dotted line). Therefore the 


algorithm used is a good representative of a BPSK signal. 
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Figure 13 - PSD of a test BPSK signal plus interference 
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Figure 14 - Theoretical and calculated PSD of a BPSK signal 


6- BLOCK S. FIL.m 


The purpose of this block is to take the FFT of the input signal 


and multiply it with the detector array filter. As stated 


3:3 


before, the photodiode array with the small space between each of 
the photodiodes can produce the filter. Placing a notch where the 
interference is present produces a narrowband frequency rejection 
via a notch filter. Also the harmful effect of notches on the other 
frequency components of the original signal must be considered. 

A plot of the filter transfer function and the test signal spectrum 
will be presented in order to visualize the effect of the filter on the 
signal. 

Finally a comparison between the filtered signal spectrum and 
the original signal spectrum is showed to evaluate a priori the 


filtering effect. 
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Figure 15 - Filter transfer function and test signal spectrum 


An example of these plots are illustrated in Figures 15 and 16. 


In Figure 15 the interference dominates the spectrum and the BPSK 
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signal appears relativelly small. Also the magnitude of the filter 
response of the excisor-detector array is shown. The overlap of the 
interference with the third notch is clearly shown. After the filtering 
process takes place (Figure 16b), the BPSK signal increases with 
respect to the interference signal, showing the effectiveness of the 


fillering process. 
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Figure 16 - Filtered and original signal spectrum comparison 
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7- BLOCK [STH 


Here the output signal is determined by an inverse transform 
of the filtered signal calculated in the previous block. In order to 
perform the inverse fast Fourier transform (IFFT) the following 


algorithm was used [Ref. 4]: 


N- 1 
s = Y, X(k +1) Wo 

Ni (20) 
The recovered signal is presented and compared with the 


original test BPSK signal without interference as shown in Figure 17. 
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Figure 17 - Original and recovered BPSK signal 


We can see that the recovered signal experiences a distortion 
due to the notches and the windowing effect of the excisor-detector 


array, in addition to the interference signal. 
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8 - BLOCK COD_B.m 


The optimum detection of a BPSK signal consists of coherent 
detection as showed in Figure 18. The gated integrator that is turned 
on at the beginning of each bit period (Tb) and integrates to the end 
of the bit. Then the integrator output is sampled and mapped into 
digital values (a binary 1 if the sample value is a positive voltage 
and a binary O if the sample value is negative). The integrator is 
then reset to a zero initial condition and this procedure is repeated 
for detection of the next bit. This optimum detector is called an 


integrate and dump detector [Ref. 7] presented in Figure 19. 


Gated 
Integrator 


cos Oc 





Figure 18 - Coherent detection 
The operation described is simulated with a FOR loop that 


computes the cumulative sum during the bit time and is dumped to 


zero at the end. To map the integrator output (Figure 20) into 
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digitals values, another loop with logic statements is performed, 
reading the middle point of each bit time. These operations generate 
a string of binary data that will be compared to the original random 


binary information as shown in Figure 21. 


Output 
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Clocking Signal 





Figure 19 - Integrate and dump detector 
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Figure 20 - Integrator output (Tb = 0.78 Msec) 
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Figure 21 - Original and recovered binary data. 


In order to detect if an error is present, the original random 
binary data and the recovered binary data are compared. This 
information is presented on the screen showing a l if no error is 


detected otherwise a -1 represents an error (Figure 22). 
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Figure 22 - Error detection 
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If in the input data the number of intervals (I) selected is 
greater than l, the program regenerates a test BPSK signal with the 
interference. frequency shifted an amount equal to the frequency 
interval (delta) introduced at the INP B.m block. This action 
simulates the movement of the excisor-detector array in order to 


evaluate the effect of this translation on the output signal. 
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An evaluation of the computer program presented in Chapter IV 


and Appendix B was performed under different conditions as shown 


in Table II. 


test agreed with what we expected or with previous experimental 


results. 


A- TEST N' 1 


For the purpose of the test, a BPSR signal vvithout interference 


was consider. 


window 9.75 MHz wide using an attenuation factor of 1 (0 dB)) was 


V. PROGRAM EVALUATION 


The main idea was to be sure that the output for each 


TABLE Il - TEST CONDITIONS 


Excisor- 
detector 


BPSK signal 
without 
interference 


BPSK signal 
with 
interference 


An excisor-detector array without notches (ideal, as a 


YES 


with notches 


Degradation of 


the original 
signal due to 
notches 


Filtering 
capability over 
the interference 
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NO 


vvithout notches 


Degradation of 
the signal 
passing through 
the window 


Degradation of 
the signal. 
Errors on the 
Data recovered 





used to represent the detector array centered at 70 MHz. The BPSK 
signal generated as a test signal had the following parameters that 


remained constant during all other tests. 


TABLE III - TEST SIGNAL PARAMETERS 


EM Amplitude BPSK signal I volt 


d Amplitude interference J volts 


Center Freq. Interference 67.4 MHz 
PRNG seed number 236578 


In the first test the central frequency of the interference was 





A 
Ai 
T 


selected as fi=0 in order to simulate a no-interference. The value of 
the bit time (Tb) equal to 0.39 psec was employed to put the minimal 
amount of the original signal on the window (main lobe and the first 
side lobes). The results of this simulation are summarized below 


(Figure 23) and match what was expected. 
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d- Recovering the 
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Figure 23 - Output Test N° 1 


As was expected the degradation of the test BPSK signal due to 
the window 9.75 MHz wide was very small and the data information 


was recovered without any error. 


B- TEST N? 2 


The second test used the same test BPSK signal as before with 
the addition of an interference. The excisor-detector was still a 
window 9.75 MHz wide with an attenuation factor of 1 (0 dB). Figure 


24 shows the results of Test N? 2. 
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a- PSD of the test 
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is 30 
db above the test 
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b- Filter signal 
and spectrum of 
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that the 
interference 
signal dominates 
the spectrum 


when a linear 


scale is used. 
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c- Output after 
passing thc test 
signal through 


the window. 


d- Recovering the 
BPSK signal and a 
comparison with 
the original test 
signal without 
interference is 


shown. 
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Figure 24 - Output Test N? 2 


As seen in Figures 24f and 24g, the interference signal produces 
a significant distortion on the recovered BPSK signal as well as on the 
random binary data. Hence a filter is desired to reduce the harmful 


interference. 
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C- TEST N* 3 


From the previous results an effective filtering process is 
needed. The actual excisor-detector array will have notches in the 
frequency domain. We want to examine the effects of these notches 
on the BPSK signal alone. The present test was performed with a test 
BPSK signal without interference passing through an excisor-detector 
with notches 60 KHz wide and with a depth of 20 dB (see Figure 


£50). 
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Figure 25 - Output Test N2 3 


No significant distortion was observed between the results of 
this filtered signal and the unfiltered BPSK signal (Test N? 1). It is 
reasonable to state that the effect of these notches (60 KHz wide and 


20 dB depth) is negligible on the broadband signal. 
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D- TEST N: 4 


Finally we consider a common situation. The test BPSK signal 
was considered with an interference present and the excisor-detector 
array had notches between detectors that will produce the filtering 
effect. A notch was aligned on the interference frequency (see Figure 


26b). 
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Figure 26 - Output Test N? 4 
During this test the effect of the narrowband notches was 


evaluated showing that this filtering technique was effective against 


the narrowband interference signal considered. 
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E- OTHER TESTS 


1- NO-FILTER TEST 


An early test was performed to check the computer program. 


The test signal was the BPSK signal without interference. No filtering 


process was 


included. We only generated and recovered the BPSK 


signal (Appendix B, p. 66). The signal was recovered without any 
distortion (Figs 27, 28 and 29). 
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Figure 27 - Original and recovered BPSK signal 
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Figure 28 - Original and recovered binary data 
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Figure 29 - Error detection output 


2 - MOVING THE EXCISOR-DETECTOR 


One of the objectives was to simulate motion of the detector 
array to study the effects of the notch passing over the interference 


signal. For this situation the program was provided with a loop that 
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simulated this movement. On the following figures the filtering 
effect can be seen. As the interference passes through the notch (60 


KHz wide and 20 dB depth), different attenuation effects take place. 
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Figure 30 - Test signal spectrum (notch at 67.35 MHz) 
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Figure 31 - Test signal spectrum (notch at 67.38 MHz) 
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Figure 32 - Test signal spectrum (notch at 67.39 MHz) 
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Figure 33 - Test signal spectrum (notch at 67.40 MHz) 
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Figure 34 - Test signal spectrum (notch at 67.41 MHz) 
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Figure 35 - Test signal spectrum (notch at 67.42 MHz) 


519 


As we can see in Figures 30 to 35, a good filtering effect is 
obtained using a relatively narrow notch with an reasonable 
attenuation. The interference signal that initially dominates the 
spectrum was progressively attenuated as the notch was moved over 


it (see Figure 32). 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This thesis has presented the basic principles, characteristics and 
potential capabilities of optical excisors. The program developed 
using MATLAB386 provides a consistent tool to aid in the study of 
acousto-optic spectrum analyzer and excision techniques. 

Processing broadband signals requires the removal of strong, 
narrowband interfering signals. Narrowband interference has been 
eliminated using a Bragg cell receiver with an excisor-detector at the 
frequency plane. Moving the excisor-detector by optical and 
mechanical means made it possible to place a notch with a given 
characteristic over the narrowband interference signal. The 
temporal characteristics of the broadband signal are recovered by 
imaging the filtered light distribution and a reference beam onto a 
photodetector array to allow coherent detection. 

All the results produced with the program during different tests 
agreed with the values obtained from available references. 

Areas for potential use of the program developed are the 
following: 

l. Evaluation of the excision effect varying the notch width and 
depth using the bit error rate as an indicator of its effectiveness. 

2. Testing different input signals and narrowband interference 


signals. 


61 


APPENDIX A - NOTATION USED IN THE PROGRAM 


The notation used consistently throughout this thesis is 


summarized for the mayor abbreviations and symbols. 


Fdb = 


Psfl 
Pssl 
PoS = 
ps» — 


Total time 

Binary random voltage (Vector) 

Signal amplitude 

Number of binary data 

Integrator Output (Vector) 

Increment (PRNG) 

Freguency interval 

Coherent detection (Vector) 

Frequency axis (Vector) 

Carrier frequency 

Interference central que nes 

Channel (Vector) 

Detector array filter signal in db (Vector) 

Detector array filter signal (Vector) 

PSD of a BPSK signal theoretical norm.(Vector) 
PSD of BPSK signal theoretical in dB (Vector) 
Recovered signal mapped into binary data (Vector) 
Number of loops in the program 

Generic counter 

Number of samples per bit time (Tb) 

Modulus (PRNG) 

Total number of elements in a vector 

Right half of a notch (Vector) 

PSD of filtered signal normalized (Vector) 

PSD of signal plus interference normalized (Vector) 
PSD of signal plus interference in dB (Vector) 

PSD of signal without interference normalized (Vec.) 


PSD of signal without interference in dB (Vector) 
Left half of a notch (Vector) 

Recovered signal in time domain (Vector) 
Random number (0 - 1) 

BPSK signal (Vector) 

Integrator and dump 

Interference signal (Vector) 

Signal (BPSK) plus interference (Vector) 

Seed for the pseudorandom number generator 
FFT of BPSK signal without interference (Vector) 
FFT of BPSK signal plus Interference (Vector) 
Signal filtered (Vector) 

Time axis (Vector) 

Generic counter 

Remainder after division (PRNG) 
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APPENDIX B - COMPUTER PROGRAM 


The following are the different blocks that using MATLAB386 as 
a master program make possible to perform the simulation and 
evaluate the results. The titles of each block in bold typeface are not 
part of the program. 


BLOCK Input initial Data (INP_B.m) 


echo on 

clc 

9o THIS IS A PROGRAM DEVELOPED IN ORDER TO EVALUATE 

% THE FILTERING CAPABILITY OF A NOTCH PRODUCED BETWEEN 
% CHANNELS IN A DETECTOR ARRAY. 


o; NOW YOU WILL BE REQUESTED TO INTRODUCE SOME PARAMETERS 


pause % Strike any key to continue 


cle 

echo on 

% THE CARRIER FREQUENCY MUST BE SELECTED BETWEEN 65 AND 75 MHz 
% WE SUGGEST 70 MHz IN ORDER TO WORK IN THE CENTER OF THE 

% DETECTOR ARRAY. 


echo off 
input Which is the Carrier Frequency (Hz)? '), 
fc=ans; % Carrier frequency 


clic 

echo on 

o THE AMPLITUDE OF THE SIGNAL MAY BE VARIED 
o WE SUGGEST TO USE 1 volt. 


echo off 

input('Select the Amplitude of the Signal = '); 

Ac=ans; % Signal amplitude 
cle 

echo on 


%. THE FREQUENCY OF THE INTERFERENCE SIGNAL MUST BE IN THE 
o FREQUENCY RANGE OF THE DETECTOR ARRAY OTHERWISE WILL BE 
% AUTOMATICALLY FILTERED (65 TO 75 MHz). 
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echo off 


input( Select the Interference Frequency (Hz) = '); 

fi=ans; % Interference Frequency 
cle 

echo on 


% IN ORDER TO RUN THE PROGRAM AND TO SIMULATE THE MOVEMENT 
% OF THE DETECTOR ARRAY YOU HAVE TO SELECT THE STEP IN 

Jo FREQUENCY THAT YOU WANT TO CONSIDER FOR EACH LOOP 

o FOR EXAMPLE: YOU WANT TO MOVE THE DETECTOR ARRAY IN 

% STEPS OF 10 KHz, INSERT THIS VALUE IN Hz (.10000 ) 


echo off 
input(‘Select the frequency interval (Hz) = ^) 
delta=ans; % Frequency Interval 


input((How many intervals do you want? '); 


I=ans; 


clc 
echo on 


% Number of loops 


% TO PROVIDE THE PROGRAM WITH ONE ESSENTIAL CHARACTERISTIC 
% IN TEST AND EVALUATION (IT MUST BE REPRODUCIBLE) IS NECESSARY % 
TO INTRODUCE A SEED FOR THE PSEUDORANDOM NUMBER GENERATOR 


(PRNG). 

echo off 

input(Input a SEED for the Random Number Generator = 5) 
seed=ans; % SEED for PRNG 

echo on 


save temp_INPUT fc fi Ac delta I 
M_BPSK % 


end 


BLOCK Main Control (M_BPSK.m) 


echo off % 

clear % 

FILTER % 

clear 

cle 

echo off 

load temp. INPUT 

for k=0:I % 

load temp_INPUT 

fc=fc+delta*k % 
% 

fi=fi+delta*k go 

pause(3) 
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Calling the main control block 


Turns program display off 
Clears memory 
Develop filter signal 


Loop for frequency steps 


Carrier frequency plus 
interval 


Interference frequency 


S. BPSR 7o Generate the input signal 


clear 7o Clears memory 

FFT B 7o Develops FFT of input 
% signal 

clear % Clears memory 

STFIL 7o Calculate the resulting 


% signal passing through the 
% detector array 


clear % Clears memory 
IFFT_B % Develops the IFFT of the 
% filtered signal 
clear % Clears memory 
COD_B % Coherent detection 
clear % Clears memory 
end % End loop 
end % End of main program 
cla % Clears program from memory 


BLOCK Detector Array Filter Signal (FILTER.m) 


load temp_S 
axis; 
cle 


% Is time to create the detector array filter signal 


goFlL-[ones(1,32768)]; % Simulates no-filter 
QS del dd: 
% left 1/2 notch 


O=[1 11111111111); 
go right 1/2 notch 
F=(Q ones(1,126) O]; 


FIL=[zeros(1,13028) 0:1:1 ones(1,137) .. 
OFFFFFFFFFFF.. 

O ones(1,137) 1:-1:0 zeros(l,.. % Detector array 
2808) 0:1:1 ones(1,137) .. % Filter signal 
OFFFFFFFFFFF.. 

O ones(1,137) 1:-1:0 zeros(1,13028)]; 

pause(3) 


cle 
o We need to create a frequency axis 


f=163.84e6*(0:16383)/32768; % Frequency axis 
clc 

% The following is a plot of the detector array 

% Filter Signal 

axis([65e6 75e6 O 1.2]); 
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plot(f,FIL(1:16384)),.. 

tile( DETECTOR ARRAY FILTER SIGNAL’)... 
xlabel FREQUENCY, H2)... 

ylabel( AMPLITUDE, volts'),.. 

pause(5) 


clg 
cle 
% Now we look at the filter magnitude response 


FIL1=FIL+eps; 
FdBz10*log10(FIL 1); 


axis([64e6 76e6 -70 10]); 

plot(f, FdB(1:16384)),.. 

title(FILTER MAGNITUDE RESPONSE’)... 
xlabel(FREQUENCY, Hz’)... 
ylabel(IH(f)l, dB’)... 

pause(5) 


clg 
cle 
save temp. F FIL f 
end 


BLOCK BPSK Signal Generator (S_BPSK.m) 


echo off 

load temp_F 

clc 

a=0.0002; % Total time 

t=0:a/32767.5:a; % Time axis 

echo on 

cle 
| % Ín order to produce a BPSK signal we need to generate 
| % a binary random voltage with values +1 or -1. 





c=16807; % Multiplier 
mo=2431-1; % Modulus 
| zi=seed; % Remainder after division 
| n=32768; % Total # of elements in the vector 
| m=64; % # of samples per bit time 
b=512; % # of binary data 
Asones(l,n), % Base vector 
p=0; 
for kzm:m:n % Loop 
g=rem(c*zi,mo); Zo Pseudorandom number 
RN=g/mo; % generator 
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Zi-g; 


d=(2*round(RN))-1; o [ +1 or -1] 
A(1,(p*m)+1:k)=A(1,(p*m)+1:k)*d; 
p=p+1; 
end % End of Loop 


% Now is time to plot the binary data 


axis([0 a -Ac-1 Ac+1]); 
plot(t,A),.. 

title(BINARY RANDOM DATA'),.. 
xlabel(' TIME, sec'),.. 

ylabel( AMPLITUDE, volts’)... 
pause 


cle 
7o Next we can form a BPSK Signal with a carrier frequency (fc) 


s=Ac*A.*cos(2*pi*fc*t); % BPSK Signal 

si-3 “Ac sin(2*pi*t^ fi), % Interference 

sd=s+si; % BPSK signal plus Interference 
cle 


% Now we are ready to plot the pulse and the BPSK signal 


axis([0 0.00002 -Ac-1 Ac+1]); 
plot(t.s(1:32768).t-AC1:32706), =. jie. 
title (BPSK TEST SIGNAL)... 
xlabel('TIME,sec')... 

ylabel( AMPLITUDE, volts'),.. 
pause 


clg 
cle 
% Plot the BPSK signal plus interference 


axis([0 0.00002 -Ac-1 Ac+1]); 

plot(t,sd),title'BPSK SIGNAL PLUS INTERFERENCE'),.. 
xlabel( TIME,sec?),.. 

ylabel('AMPLITUDE ),.. 

pause 


clg 
clc 


save temp S mts Ac fcasdAnb 


save temp. SR m t AC A n b a 
end 
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BLOCK FFT of BPSK Signal (FFT_B.m) 


echo on 
load temp_S 
load temp_F 


cle 
9o Now ve get the Fourier transform of the input signal (BPSK) 


S-fft(s); % FFT BPSK Signal 
SD=fft(sd); % FFT BPSK Signal + Interference 
cle 


% In order to plot the power spectral density of the signal we need to 
% make some calculations 


Pss=SD.*conj(SD); 


Pssl-Pss/max(Pss); % PSD Signal + Interference 
PSS=10*log10(Pss1); % PSD in dB 

Psl=S.*conj(S); 

Ps2=Ps1/max(Ps1); % PSD Signal without interference 
PSS 1=10*log10(Ps2); % PSD in dB 

cle 


% Plot the power spectral density of the BPSK signal 
7o plus interference 


axis([fc-5e6 fc+5e6 -90 0)); 

plot(f,PSS(1:16384)),.. 

title(PSD OF BPSK SIGNAL PLUS INTERFERENCE ).,.. 
xlabel FREQUENCY, Hz)... 

ylabel POWER, dB)... 

pause 

hold on 


clc 
% For comparison we calculate the theoretical PSD of a BPSK signal 


R=(1/((a/32768)*m)); 
L=(pi*(f-fc)/R)+eps; 
K=(Ac‘2)/(4*R); 
Ne(sin(L)./L).^2; 
PSD=K*X; 
G=PSD./max(PSD); 
G1-10*log10(G); 


clc 
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% Now we can plot the theoretical PSD of the BPSK 


plot(f,G1(1:16384)),.. 
meta GPSDB,pause 
hold off 


clg 

cele 

% For comparison purpose we plot the BPSK signal calculated 
% without interference and the theoretical one. 


axis([fc-5e6 fc+5e6 -90 0)); 
plot(f,PSS1(1:16384)),.. 

title(BPSK SIGNAL COMPARISON’)... 

xlabel (FREQUENCY, Hz)... 

ylabel POWER, dB)... 

text (0135008505 —— 08 THEORETICAL (dotted line)’,.. 
text(0.13,0.80,' CALCULATED (solid line)',.. 
pause 

hold on 

plot(f,G1(1:16384)),.. 

pause 

hold off 


clg 
cle 


save temp. FFT SD f FIL Pssl S fc 
end 
BLOCR Signal Filtering on the Detector Array (S FIL.m) 
echo on 
load temp. FFT 


clc 
% Now is time to pass the BPSK signal through the detector array 


SF=SD.*FIL; % Filtering the Signal 
Psf=SF.*conj(SF); 

Psf1=Psf/max(Psf); % PSD Filtered Signal 
cle 


% First of all we compare the signal spectrum with the filter due 
7o to the detector array 


axis([fc-6e6 fc+6e6 0 1.2]); 


plot(f,Pssl(1:16384),f,FIL(1:16384),':'),.. 
title(FILTER SIGNAL & BPSK SPECTRUM ^)... 
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xlabel (FREQUENCY, Hz)),.. 
ylabel( AMPLITUDE, volts'),.. 
meta GSFIL,pause 


clg 
cle 
% Now we can plot the signal filtered compared to the original one 


axis([fc-6e6 fc+6e6 O 1.2]); 
plot(f,Psf1(1:16384),f,Pss1(1:16384),':'),.. 
tile(BPSK SPECTRUM COMPARISON)... 
xlabel( FREQUENCY, Hz)),.. 

ylabel( AMPLITUDE, volts'),.. 

pause 


clg 
elc 


save temp. SF SF f 
delete temp, FFT.mat 
end 


BLOCK IFFT Signal Filtered (IFFT B.m) 


echo on 
load temp. SF 
load temp S 


clc 
% Now we return to the time domain taking the inverse 
% Fourier transform of the filtered signal 


rszifft(SF); % Inverse Fourier transform 
t=1'; % Arrange time axis 
clc 


% We can plot the Recovered Signal in the time domain 


axis([0 0.00002 -Ac-1 Ac+1]) 
plot(t,rs) title (RECOVERED BPSK SIGNAL )),.. 
xlabel( TIME, sec?),.. 

ylabel AMPLITUDE, volts'),.. 

pause 


| 
i 


clg 
clc 
7o Let's compare with the original signal plotting it 


71 





subplot(211),plot(t,s); title (ORIGINAL BPSK’)... 

xlabel( TIME, sec'),ylabel (AMPLITUDE'),.. 
subplot(212),plot(t,rs);title/'RECOVERED BPSK SIGNAL )... 
xlabel(oTIME, sec )ylabel AMPLITUDE),.. 

pause 


clg 

Cre 

save temp_R rs ts m Ac fc 
delete temp_SF.mat 

delete temp_S.mat 

end 


BLOCK Coherent Detector (COD_B.m) 


echo on 

load temp_R 

load temp_SR 

clc 

% Now is time to make a coherent detection. The optimum 
% detection consist of product detection with an integrator 
% that is turned on at the beginning of each bit period 

% and integrates to the end of the period. 


DET—zrs.*cos(2*pi*fc*t); % Coherent detection 
B=ones(1,n); % Base vector 
for u=0:m:n-m; % Loop 
sa=cumsum(DET(1,(u+1:u+m))); % Integrator 
B(1,(u+1:u+m)=B(1,(u+1:u+m)).*sa; 
sa=0; % Dump 
end % End loop 
P=-1; 
Hzones(l,n); 
for v=0:m:n-m; % Loop 
dig=B(1,v+m/2); % Reading the middle 
point 
if dig>0 
H(1,(v+1:v+m))=H(1,(v+1:v+m)).*1; % Mapped into a digital 
value 
elseif dig<O 
H(1,(v+1:v+m))=H(1,(v+1:v+m)).*P; % Mapped into a digital 
value 
end % End conditional 
Statement 
end % End Loop 


UE 


E-A.*H; % Checking for error 


cle 
% We look at the output of the integrator 


axis({0 a -m-1 m+1)); 
plot(t,B),.. 

title INTEGRATOR OUTPUT’)... 
xlabel(TIME, sec'),.. 

ylabel( AMPLITUDE, volts'),.. 
pause 


clg 

clc 

% Comparison of the Random Binary Data and the data 
% recovered 


axis([Ü a -Ac-2 Ac42]); 

subplot(211),plot(t,A); title (ORIGINAL BINARY DATA)... 
xlabel(' TIME, sec'),ylabel( AMPLITUDE'),.. 
subplot(212),plot(t,H);title(RECOVERED BINARY DATA"... 
xlabel('TIME, sec'),ylabel( AMPLITUDE'),.. 

pause 


clg 
elc 


76 In order to determine if an error is present we compare 
% the data. If we get all 1, there is "no error"; otherwise a -1 
% represents an "error". 


axis([0 a -2 2J); 

plot(t,E),title( ERROR DETECTION")... 
text(0.5,0.75, NO ERROR','sc”),.. 
text(0.5,0.25, ERROR','sc'),.. 

xlabel(' TIME, sec'),.. 

ylabel( AMPLITUDE, volts')... 
pause 


clg 
clc 


delete temp S.mat 


delete temp R.mat 
end 
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APPENDIX C - INFORMATION ABOUT NOTCHES 


In order to help the user in the selection of the interference 
frequency and evaluate the effect in a particular notch, here is a 
table showing the position of each notch depending how wide it is in 


the frequency domain. 


Notch Notch Notch 10 KHz Notch 20 KHz 

“ Center Start freq. End freq. Start freq. End freq. 
l 65900000 65895000 65905000 65890000 65910000 
2 66650000 66645000 66655000 66640000 66660000 
3 67400000 67395000 67405000 67390000 67410000 
4 68150000 68145000 68155000 68140000 68160000 
5 68900000 68895000 68905000 68890000 68910000 
6 69650000 69645000 69655000 69640000 69660000 
7 70400000 70395000 70405000 70390000 70410000 
8 71150000 71145000 71155000 71140000 71160000 
9 71900000 71895000 71905000 71890000 71910000 
10 72650000 72645000 72655000 72640000 72660000 
1] 73400000 73395000 73405000 73390000 73410000 
WA 74150000 74145000 74155000 74140000 74160000 
Notch Notch Notch 30 KHz Notch 40 KHz 

a Center Start freq. End freq. Start freq. End freq. 
l 65900000 65885000 65915000 65880000 65920000 
2 66650000 66635000 66665000 66630000 66670000 
3 67400000 67385000 67415000 67380000 67420000 
4 68150000 68135000 68165000 68130000 68170000 
5 68900000 68885000 68915000 68880000 68920000 
6 69650000 69635000 69665000 69630000 69670000 
7 70400000 70385000 70415000 70380000 70420000 
8 71150000 71135000 71165000 71130000 71170000 
9 71900000 71885000 71915000 71880000 71920000 
10 72650000 72635000 72665000 72630000 72670000 
11 73400000 73385000 73415000 73380000 73420000 
12 74150000 74135000 74165000 74130000 74170000 


Notch 
H 


OO JO A B UY N = 


mi pė pė 
N — © 


Notch 


NO 00 IJ ON Mr RW IN = 


Notch 


NO 00 IJ Gh Ln E LM t9 — 


ph ph ked 
N — © 


Notch 
Center 


65900000 
66650000 
67400000 
68150000 
68900000 
69650000 
70400000 
71150000 
71900000 
72650000 
73400000 
74150000 


Notch 
Center 


65900000 
66650000 
67400000 
68150000 
68900000 
69650000 
70400000 
71150000 
71900000 
72650000 
73400000 
74150000 


Notch 
Center 


65900000 
66650000 
67400000 
68150000 
68900000 
69650000 
70400000 
71150000 
71900000 
72650000 
73400000 
74150000 


Notch 50 KHz 


Stan freq. 


65875000 
66625000 
67375000 
68125000 
68875000 
69625000 
70375000 
71125000 
71875000 
72625000 
73375000 
74125000 


End freg. 


65925000 
66675000 
67425000 
68175000 
68925000 
69675000 
70425000 
71175000 
71925000 
72675000 
73425000 
74175000 


Notch 70 KHz 


Start freg. 


65865000 
66615000 
67365000 
68115000 
68865000 
69615000 
70365000 
71115000 
71865000 
72615000 
73365000 
74115000 


End freg. 


65935000 
66685000 
67435000 
68185000 
68935000 
69685000 
70435000 
71185000 
71935000 
72685000 
73435000 
74185000 


Notch 90 KHz 


Start freg. 


65855000 
66605000 
67355000 
68 105000 
68855000 
69605000 
70385000 
71105000 
71855000 
72605000 
73355000 
74105000 


End freg. 


65945000 
66695000 
67445000 
68195000 
68945000 
69695000 
70445000 
71195000 
71945000 
72695000 
73445000 
74195000 


Notch 60 KHz 


Start freg. 


65870000 
66620000 
67370000 
68120000 
68870000 
69620000 
70370000 
71120000 
71870000 
72620000 
73370000 
74120000 


End freg. 


65930000 
66680000 
67430000 
68 180000 
68930000 
69680000 
70430000 
71180000 
71930000 
72680000 
73430000 
74180000 


Notch 80 KHz 


Start freq. 


65860000 
66610000 
67360000 
68110000 
68860000 
69610000 
70360000 
71110000 
71860000 
72610000 
73360000 
74110000 


End freq. 


65940000 
66690000 
67440000 
68190000 
68940000 
69690000 
70440000 
71190000 
71940000 
72690000 
73440000 
74190000 


Notch 100 KHz 


Start freq. 


65850000 
66600000 
67350000 
68100000 
68850000 
69600000 
70350000 
71100000 
71850000 
72600000 
73350000 
74100000 


End freq. 


65950000 
66700000 
67450000 
68200000 
68950000 
69700000 
70450000 
71200000 
71950000 
72700000 
73450000 
74200000 


Notch 
# 


RO QO d C Ln fh WW — 


peu beet jm 
N © 


Notch 
Center 


65900000 
66650000 
67400000 
68150000 
68900000 
69650000 
70400000 
71150000 
71900000 
72650000 
73400000 
74150000 


Notch 110 KHz 
Start freq. 


65845000 
66595000 
67345000 
68095000 
68845000 
69595000 
70345000 
71095000 
71845000 
72595000 
73345000 
74095000 


End freq. 


65955000 
66705000 
67455000 
68205000 
68955000 
69705000 
70455000 
71205000 
71955000 
72705000 
73455000 
74205000 


Notch 120 KHz 


Stan freq. 


65840000 
66590000 
67340000 
68090000 
68840000 
69590000 
70340000 
71090000 
71840000 
72590000 
73340000 
74090000 


End freq, 


65960000 
66710000 
67460000 
68210000 
68960000 
69710000 
70460000 
71210000 
71960000 
72710000 
73460000 
74210000 


To put the correct value in the notch filter that represent a given 
attenuation in dB, a table with different values of attenuation from O 


to -45 dB is presented. 


ATTENUATION TABLE 


Attenuation Number Attenuation Number 


dB 


1.00000 
0.79432 
0.63096 
0.50119 
0.39811 


0.31623 
0.25119 
0.19953 
0.15849 
0.12589 
0.10000 
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dB 


0.07943 
0.06309 
0.05119 
0.03981 
0.03162 
0.01000 
0.00316 
0.00100 
0.00032 
0.00010 
0.00003 





APPENDIX D - DEFINITIONS 


A group of useful definitions related to this work are presented. 


Bragg cell : The heart of all acousto-optic systems. A thin slab 
of a transparent crystal (such as tellurium oxide, lithium niobite, 
or gallium phosphide) with one or more piezoelectric transducers 
attached. When a microwave signal excites the transducer, it 
creates an ultrasound wave into the crystal. This wave changes 
the optical index of refraction at the microwave frequency [Ref. 


DI. 


Excisor-detector: An array of photodetectors which is inserted 
into the frequency plane of an optical spectrum analyzer. The 
small spaces between detectors are used as spatial light 


modulators, that produce notches in the frequency response. 


Model : Mathematical representation of a system or part of a 
system and/or some portion of its environment. The model is a 
set of equations which may be solved manually or by means of a 
computer to determine the behavior of the system under a 
specific set of conditions or during a series of sets of conditions 
Models permit easy modification of the system parameters to 
determine the effects of such changes on system performance 


(Ref. 10). 


CRE 


Simulation : Simulations include techniques for representing a 
system Or part of a system and/or some portion of its 
environment (including human interfaces with hardware) 
through the use of sustitutes for all or part of the real hardware, 
personnel; and environment. Thus simulation covers a more 
extensive field than modeling, although most models are also 


simulations [Ref. 10]. 
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